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This  report  describes  the  RENDER  calculatlonal  procedure  that  was  de''el- 
■oped  for  use  in  evaluating  the  spatial-  and  time-dependent  energy  deposition 
' in  air  produced  by  anisotropic  nuclear  sources.  The  RENDER  procedure  utilizes 
nuclear  energy  deposition  data  generated  for  conical  sources  that  were  origin- 
ally computed  for  line-beam  sources  through  use  of  the  ZAP,  ZAPN,  ZAPGAM  and 
) DEPO  codes.  The  energy  deposition  data  for  air  that  are  input  to  RENDER  were 
^ generated  for  neutron  and  gamma-ray  sources  and  for  secondary  gamma  rays  re- 
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suiting  from  neutron  capture  and  inelastic  scattering  in  air.  The  energy 
deposition  data  for  line-beam  sources  and  for  point  isotropic  sources  were 
found  to  compare  favorably  with  similar  data  reported  in  the  literature. 

The  RENDER  procedure  was  run  utilizing  energy  deposition  data  from  the  conical 
source-data  base  for  a 9-to-lO-MeV  gamma-ray  source  and  the  results  of  the 
convolution  over  source  emission  direction  and  time  were  found  to  be  in  good 
agreement  with  the  input  data,  indicating  that  the  RENDER  procedure  performs 
the  tlme-and-angle  convolution  correctly^ 

I 

Volumes  II  through  IV  of  this  report  present  tabulated  data  on  the  time  depen-j 
dent  energy  deposition  in  air  versus  range  for  neutron  and  gamma-ray  point 
Isotropic  sources  and  for  secondary  gamma  rays  generated  by  point  Isotropic 
I neutron  sources.  Also  given  in  Vol.  V are  curve  fit  coefficients  for  use  in 
computing  the  energy  deposition  in  air  versus  distance  and  source  emission 
angle  for  line  beam  sources  of  neutrons  and  gamma  rays.  Coefficient  data  are 
also  given  for  secondary  gamma-ray  energy  deposition  by  line  beam  neutron 
sources.  | 
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I . INTRODUCTION 


The  studies  described  In  this  report  were  performed  for  the  purpose 
of  providing  a calculatlonal  procedure  for  predicting  the  time  and  spatial 
dependency  of  the  nuclear  energy  deposition  (neutron,  gamma  ray,  and  neutron- 
produced  secondary  gamma  rays)  In  an  Infinite  homogeneous  medium  of  air  for 
anisotropic  nuclear  sources.  There  were  no  data  available  In  the  literature 
at  the  time  this  study  was  started  that  would  provide  the  needed  source 
energy,  time-  and  spatial-dependent  nuclear-energy  deposition  data  for  use 
In  the  calculatlonal  procedure. 

The  time-  and  spatial-dependent  nuclear-energy  deposition  calcula- 
tions were  performed  as  a function  of  source  energy  first  for  line-beam 
neutron  and  gamma-ray  sources  and  the  resulting  data  were  then  converted 
to  give  energy  deposition  data  for  conical  sources.  A computer  procedure, 
RENDER,  was  developed  to  use  the  energy  deposition  data  for  conical  sources 
to  calculate  the  source  energy,  time-  and  spatial-dependent  energy  depo- 
sitions resulting  from  anisotropic  nuclear  sources. 

The  atmosphere  was  assumed  to  be  a homogeneous  mixture  of  nitrogen 
and  oxygen  with  a density  of  1.225  grams  per  liter.  The  nuclear  cross 
sections  used  In  the  calcualtlons  are  described  In  Section  II.  The  nuclear 
energy  deposition  calculations  were  performed  using  the  ZAP,  ZAPN  and  ZAPGAM 
Monte  Carlo  procedures  and  the  DEPO  code  (Refs.  1,  2,  3,  and  4).  A discus- 
sion of  the  calculatlonal  methods  and  the  results  obtained  for  line-beam 
gamma-ray  and  neutron  sources  and  the  secondary  gamma-ray  sources  resulting 
from  gamma-ray  production  by  neutron  Inelastic  scattering  and  n,Y  reactions 
are  also  given  In  Section  II. 

The  time- Independent  energy-deposition  results  obtained  with  these 
codes  were  compared  with  published  data  to  assess  the  accuracy  of  the  newly- 
calculated  data.  The  results  of  these  comparisons  are  given  In  Section  II. 

The  collision  data  generated  for  line-beam  sources  with  use  of  the 
ZAP,  ZAPN  and  ZAPGAM  Monte  Carlo  procedures  were  stored  on  history  tapes. 
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2 

These  history  tapes  were  then  used  as  Input  to  a modified  version  of  the 
DEPO  procedure  to  generate  time-  and  spatial-dependent  energy-deposition 
data  for  conical  neutron  and  gamma-ray  sources.  The  calculatlonal  methods 
used  are  discussed  in  Section  III. 

A computer  procedure,  designated  as  RENDER,  was  written  to  convolute 
the  time-  and  spatial-dependent  energy-deposition  data  for  conical  sources 
with  time-  and  angle-dependent  source  data  for  anisotropic  nuclear  sources. 
The  RENDER  procedure  is  described  in  Section  IV.  Some  results  of  RENDER 
calculations  that  were  run  to  check  out  the  accuracy  of  the  calculatlonal 
methods  used  in  the  RENDER  code  are  also  described  in  Section  IV. 

Of  special  interest  to  laboratories  working  on  problems  related 
to  the  exposure  of  personnel  and  equipment  to  nuclear  radiation  is  1)  the 
time-dependent  energy-deposition  data  for  air  vs  range  for  point-isotrop- 
ic sources;  2)  the  time-integrated  energy-deposition  data  for  air  vs  range 
for  point-isotropic  sources;  and  3)  the  time-integrated  energy-deposition 
data  for  air  vs  range  and  source  polar  angle  for  line-beam  sources. 

Tables  containing  these  data  for  neutron,  gamma-ray  and  secondary  gamma- 
ray  sources  for  each  of  the  neutron  and  gamma-ray  source  energy  intervals 
ccnsldered  in  this  study  are  given  in  Vols.  II  through  V of  this  report. 
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II.  ENERGY  DEPOSITION  CALCULATIONS  FOR  LINE-BEAM  SOURCES 

For  the  ganm^-ray  transport  calculations,  the  atmosphere  was 
considered  to  be  a homogeneous  mixture  of  nitrogen,  oxygen,  helium  and  argon 
with  number  densities  of  3.978X10^^,  1.067X10^^,  1.335X10^^  and  2.380X10^^, 
respectively.  The  density  of  the  atmosphere  was  assumed  to  be  1.225  grams/ 
liter. 

The  cross  section  data  for  air  used  in  the  ZAP  and  ZAPGAM  calcula- 
tions for  the  total  cross  section,  the  photoelectric  cross  section,  the 
pair  production  cross  section, the  coherent  scattering  cross  section  and 
the  Incoherent  scattering  cross  section  as  a function  of  energy  were  taken 
from  Ref.  5 for  67  energies  between  0.1  keV  and  100  MeV.  The  coherent 
and  incoherent  scattering  structure  factors  given  in  Ref.  5 for  A5  values 

of  the  momentum  transfer  between  0.0  and  4.85X10  (M  c units)  were  also 

e 

used  in  the  gamma-ray  scattering  calculations. 

Neutron  cross  section  data  needed  for  Input  to  ZAPN  were  generated 
from  ENDF/B  -III  files  (Ref.  6)  for  nitrogen  and  oxygen  with  use  of  the 
SAM-X  and  CBAND  procedures  (Ref.  7).  The  ENDF/B  files  were  first  employed 
a’  Input  to  SAJl-X  which  extracted  basic  cross  section  data  from  the  ENDF/B- 
III  files  for  nitrogen  and  oxygen  and  organized  the  data  into  the  formats 
required  by  CBAND.  The  output  from  SAM-X  is  a binary  tape  (element  data 
tape,  ETD)  which  was  then  input  to  CBAND  which  organized  the  cross  section 
data  into  energy  bands.  The  output  of  CBAND  is  a magnetic  tape  which  is 
referred  to  as  an  organized  data  tape  (ODT) . It  is  the  organized  data  tape 
that  is  used  as  input  to  ZAPN. 

For  the  neutron  transport  problems  it  was  assumed  that  air  was  78% 
nitrogen  and  22%  oxygen  by  volume  and  the  air  density  was  taken  to  be 
1.225  grams/liter. 

In  ZAPN  the  inelastic  scattering  cross  section  is  considered  to  be 


the  sum  of  all  interactions  that  produce  an  x-partlcle  (alpha,  protons, 
gammas,  etc.)  in  addition  to  inelastically  scattering  the  incident  neutron. 
For  nitrogen,  the  inelastic  scattering  cross  section  at  energies  above 
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8 MeV  is  the  sum  of  the  (n,n'p),  (n,n'’y')  and  (n,n'a)  cross  sections  and 
thus  not  all  of  the  energy  required  to  excite  the  nucleous  to  a level  or 
into  the  continuum  (the  Q value)  goes  into  the  production  of  gamma  rays. 
ZAPN  was  modified  so  that  when  an  inelastic  collision  occurred  with  a 
nitrogen  atom,  a random  number  was  used  to  determine  if  the  Q value  was 
to  be  given  off  in  the  form  of  gamma  rays  or  if  charged-particle  emission 


occurs.  If  O , (E)  is  the  inelastic  scattering  cross  section  for  the 
n.n'Y 

"x"  particles  being  gamma  rays  and  O , (E)  is  the  total  Inelastic  scat- 

n,n 

tering  cross  section,  then 


a .^(E) 

P(E)  = — Sj-S— ^ — 

’ a ,(E) 
n,n 


is  the  fraction  of  the  inelastic  collisions  for  which  gamma  rays  are 
emitted  that  carry  an  energy  equal  to  the  Q value  and  (1  - P(E))  is  the 
fraction  of  the  inelastic  collisions  for  which  an  energy  equal  to  the 
Q value  is  carried  off  by  charged  particles,  such  as  alpha  particles  and 
protons.  When  an  inelastic  scattering  event  was  selected  to  be  a non- 
gamma-ray-producing  event,  then  an  energy  equal  to  the  Q value  was  added 
to  the  counter  denoting  the  energy  deposited  by  absorption  processes. 

For  nitrogen,  P(E)  was  calculated  from  ENDF/B-III  cross  sections.  Table 
1 lists  these  values  for  neutron  energies  between  8 and  16  MeV.  Linear 
interpolation  was  used  on  the  data  in  Table  I to  obtain  values  of  P(E) 
at  energies  between  the  tabulated  energies. 

The  neutron  history  tape  generated  with  the  ZAPN  procedure  was 
used  in  the  SGSORC  procedure  (Ref.  8)  to  generate  secondary  gamma-ray 
photon  sources.  The  number  of  gamma  rays  produced  per  neutron  collision 
as  a function  of  neutron  energy  and  the  energy  distribution  of  the  emitted 
gamma  rays  as  a function  of  the  incident  neutron  energy  were  obtained 
from  the  ENDF/B-IIl  files  for  use  in  the  SGSORC  procedure. 


TABLE  I.  FRACTION  OF  THE  NONELASTIC  COLLISIONS  WITH 
NITROGEN  THAT  RESULT  IN  GAMMA-RAY  EMISSION 


NEUTRON 

ENERGY (MeV) 

P(E)  - 

on,n 

8.00 

1.000 

8.50 

0.974 

9.00 

0.900 

9.50 

0.816 

9.75 

0.764 

10.00 

0.708 

10.25 

0.678 

10.05 

0.655 

10.75 

0.629 

11.00 

0.604 

11.25 

0.571 

11.50 

0.543 

11.75 

0.517 

12.00 

0.496 

12.25 

0.477 

12.50 

0.456 

13.00 

0.430 

13.50 

0.414 

14.00 

0.403 

14.50 

0.395 

15.00 

0.387 

15.50 

0.378 

16.00 

0.360 
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For  the  ZAPGAM  calculation  of  the  energy  deposition  by  secondary 
gamma  rays,  the  same  gamma-ray  cross-section  data  used  in  the  ZAP  calcula- 
tions were  also  used  in  ZAPGAM. 


The  geometry  used  in  the  energy  deposition  calculations  for  line- 
beam  sources  is  shown  in  Fig.  1.  The  source  radiation  (neutron  or  gamma 


ray)  is  emitted  from  the  source  point  in  the  direction  6=0*.  The  energy 
deposition  was  determined  for  regions  formed  by  the  Intersection  of  the 
conical  surfaces  generated  by  the  revolution  of  the  radii  ^ and  R^ 


inclined  at  polar  angles  Qj  > respectively,  about  the  source 

axis  and  the  spherical  surfaces  formed  by  the  revolution  of  the  radii 


R^_^  and  R^  about  the  source  point.  The  radii  and  polar  angles  used  to 
define  the  energy  deposition  regions  out  to  1500  meters  are  listed  in 


Table  II. 


Because  of  the  fact  that  gamma  rays  travel  at  the  speed  of  light 
whereas  neutrons  travel  at  a speed  determined  by  their  energies,  different 
time  intervals  were  used  to  record  the  time-dependent  energy  deposition 
data  for  gamma  rays,  neutrons  and  secondary  gamma  rays.  Table  III  lists 
the  time-interval  bounds  used  for  each  of  the  three  different  radiation 


sources. 


The  source  energy  Intervals  used  for  the  neutron  and  gamma-ray 

line-beam  sources  are  listed  In  Table  IV.  For  the  line-beam  gamma-ray 

problems,  10,000  histories  were  run  for  each  source  energy  interval  and 

histories  were  terminated  either  by  a minimum  weight  cutoff  of  1.0X10  ^ 

or  when  the  number  of  collisions  per  history  exceeded  30.  For  the  line- 

beam  neutron  problems,  2000  histories  were  run  for  each  source  energy 

interval.  Neutron  histories  were  terminated  when  the  number  of  collisions 

-3 

per  history  exceeded  200,  when  the  energy  reduced  below  1.01X10  keV  or 

when  the  transport  time  exceeded  1.0  second.  The  secondary  gamma-ray 

histories  were  terminated  when  the  number  of  collisions  exceeded  30,  when 

-3 

the  particle  weight  reduced  by  a factor  of  1X10  or  when  the  gamma-ray 
energy  had  reduced  below  0.2  keV.  Secondary  gamma-ray  sources  were  gener- 
ated by  SGSORC  using  10  per  cent  of  the  collision  data  from  the  first  1000 
of  the  neutron  histories  stored  on  the  ZAPN-generated  collision  tapes. 


% 
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CONICAL  SOURCE  GEOMETRY 


ENERGY  DEPOSITION  GEOMETRY 


Fig.  1.  Energy  Deposition  Geometry  for  Line  Beam  and  Conical  Sources 
and  Source  Geometry  for  Conical  Sources 
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TABLE  II.  RADIAL  AND  ANGLE  INTERVALS  USED  TO  DEFINE 
DEPOSITION  GEOMETRY 


Radial  Interval  Polar  Angle  Interval 

^m) (degrees) 


0 

- 

10 

0 

- 

1 

10 

- 

20 

1 

- 

2 

20 

- 

40 

2 

- 

3 

40 

- 

60 

3 

- 

5 

60 

- 

80 

5 

- 

7 

80 

- 

100 

7 

- 

10 

100 

- 

125 

10 

- 

15 

125 

- 

150 

15 

- 

20 

150 

- 

175 

20 

- 

30 

175 

- 

200 

30 

- 

40 

200 

- 

250 

40 

- 

50 

250 

- 

300 

50 

- 

60 

300 

- 

350 

60 

- 

70 

350 

- 

400 

70 

- 

80 

400 

- 

450 

80 

- 

90 

450 

- 

500 

90 

- 

100 

500 

- 

600 

100 

- 

120 

600 

- 

800 

120 

- 

140 

800 

- 

1000 

140 

- 

160 

1000 

- 

1200 

160 

- 

180 

1200 

- 

1500 

» 
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TABLE  III.  BOUNDS  OF  DEPOSITION  TIME  INTERVALS  USED  IN  NEUTRON, 
SECONDARY  GAMMA-RAY  AND  PRIMARY  GAMMA-RAY  PROBLEMS 


DEPOSITION  TIME  INTERVAL  (usec) 


Neutrons  and  Secondary 

Gamma-Ray  Problems Primary  Gamma  Problems 


0 

0.2 

0 

1 

O 

0.2 

- 

0.3 

10-^^ 

- 

0.02 

0.3 

- 

0.5 

0.02 

- 

0.05 

0.5 

- 

1.0 

0.05 

- 

0.07 

1 

- 

2 

0.07 

- 

0.10 

2 

- 

3 

0.10 

- 

0.15 

3 

- 

5 

0.15 

- 

0.20 

5 

- 

7 

0.20 

- 

0.30 

7 

- 

10 

0.30 

- 

0.50 

10 

- 

20 

0.50 

- 

0.70 

20 

- 

35 

0.70 

- 

1.00 

35 

- 

50 

1.00 

- 

1.50 

50 

- 

70 

1.50 

- 

2.00 

70 

- 

100 

2.00 

- 

3.00 

100 

- 

200 

3.00 

- 

5.00 

200 

- 

400 

5.00 

- 

7.00 

400 

- 

700 

7.00 

- 

10.00 

700 

- 

1000 

1000 

- 

2000 

2000 

- 

3500 

3500 

- 

15000 

15000 

- 

100,000 

J 
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TABLE  IV.  SOURCE  ENERCY  INTERVAL  BOUNDS  USED  FOR  NEUTRON  AND 
PRIMARY  CAMMA-RAY  PROBLEMS 


Neutron  Source  Primary  Gamma-Ray 

Energy  Intervals  (MeV)  Energy  Intervals  (MeV) 


0.001 

- 

0.00335 

9.0 

- 

10.0 

0.00335 

- 

0.0912 

8.0 

- 

9.0 

0.0912 

- 

0.0248 

7.0 

- 

8.0 

0.0248 

- 

0.0676 

6.0 

- 

7.0 

0.0676 

- 

0.184 

5.0 

- 

6.0 

0.184 

- 

0.303 

4.0 

- 

5.0 

0.303 

- 

0.  50 

3.0 

- 

4.0 

0.50 

- 

0.823 

2.0 

- 

3.0 

0.823 

- 

1 .353 

1.0 

- 

2.0 

1.353 

- 

1.738 

0.5 

- 

1.0 

OO 

- 

2.232 

0.1 

- 

0.5 

2.232 

- 

2.865 

0.01 

- 

0.1 

2.865 

- 

3.680 

3.680 

- 

6.070 

6.070 

- 

7.  790 

7.790 

- 

10.0 

10.0 

- 

12.0 

12.0 

- 

13.5 

13.5 

- 

15.0 

The  line-beam  source  problems  were  run  using  the  computer  code 
outlined  in  Fig.  2.  The  ZAPXG  procedure  was  used  to  generate  collision 
data  for  line-beam  gamma-ray  sources.  The  ZAPNOl  procedure  was  used  to 
generate  collision  data  for  line-beam  neutron  sources.  The  collision 
data  generated  by  the  line-beam  neutron  and  gamma-ray  problems  were 
stored  on  magenetic  tapes  (history  tapes).  The  collision  data  on  the 
neutron  history  tapes  were  rearranged  with  use  of  the  RETAPE  code  so 
that  the  resulting  neutron  history  tapes  could  be  read  by  both  the 
SGSROl  procedure  and  the  DEPO  procedure.  The  SGSROl  procedure  was  used 
with  the  neutron  history  tapes  to  generate  secondary  gamma-ray  source 
data  for  input  to  the  ZAPGAM  procedure  which  was  then  used  to  generate 
secondary  gamma-ray  history  tapes.  The  DEPO  procedure  was  used  to  read 
the  history  tapes  generated  by  each  of  the  ZAPXG,  ZAPNOl  and  ZAPGAM  runs 
and  to  calculate  the  energy  deposition  occurring  within  the  energy 
deposition  volumes  and  deposition  time  intervals  given  in  Tables  11  and 
III.  The  output  of  the  DEPO  runs  were  obtained  in  both  printed  and 
computer  tape  formats.  For  each  ZAP,  ZAPN  and  ZAPGAM  problem  run,  the 
DEPO  printed  output  gave  3 pages  of  printout  (energy  deposition  vs  radial 
interval  and  polar  angle  interval)  for  each  source  time  interval  and  each 
source  energy  Interval  and  a final  3 pages  for  the  time- Integrated  data 
for  each  source  energy  interval.  This  resulted  in  648  pages  of  output 
for  the  primary  gamma-ray  problems,  1197  pages  of  output  for  the  neutron 
problems  and  1311  pages  output  for  the  secondary  gamma-ray  problems. 

The  printed  output  from  the  DEPO  procedure  and  the  Monte  Carlo 
runs  were  examined  to  Insure  that  no  energy  was  lost  in  the  Monte  Carlo 
calculations.  Table  V lists  the  fraction  of  the  primary  gamma-ray  source 
energy  that  was  deposited  in  the  atmosphere  within  1500  meters  slant  range 
by  gamma-ray  scattering  and  absorption  Interactions.  The  difference  between 
1 and  the  fraction  of  the  source  energy  deposited  In  the  atmosphere  is  the 
fraction  of  the  source  energy  that  was  lost  due  to  the  use  of  Russian 
Roulette  in  terminating  histories  by  minimum  weight  cutoff,  by  photons 
that  underwent  more  than  30  collisions  per  history  and  by  photons  that 


/ 

X- Sect Ions  ^ 

ZAPXG 

Ref.  5 ^ 

Ref  • X 

History 

Tape 


X- Sect  Ions 
Ref.  5 


ZAPGAM 
Ref.  3 


History 

Tape 


SEC.  y'sj 


PRIMARY 

GAMMA-RAYS 


Deposition 

Tape 


NEUTRONS 

Deposition 

Tape 


SECONDARY 

GAMMA-RAYS 


Deposition 

Tape 


Fig.  2.  Calculation  Flow  for  Energy  Deposition  Calculations 


TABLE  V,  ENERGY  DEPOSITED  FOR  SLANT  RANGES 
LESS  THAN  1500  METERS  AND  AVERAGE  NUMBER  OF 
COLLISIONS  PER  HISTORY  FOR  GAMMA-RAY  SOURCES 


Gamma- Ray 

Source  Energy 
Interval  (MeV) 

Fraction  of  Source 
Energy  Deposited 
by  Scattering 
and  Absorption 
(MeV) 

Average  Numbe: 
of  Collisions 
per  History 

0.01-0. 1 

0.9697 

13.0 

0 

1 

o 

0.9916 

22.6 

0. 5-1.0 

0.9971 

24.4 

1. 0-2.0 

0.9965 

25.1 

2. 0-3.0 

0.9928 

25.4 

3. 0-4.0 

0.9890 

25.5 

4. 0-5.0 

0.9846 

25.6 

5. 0-6.0 

0.9788 

25.7 

6. 0-7.0 

0.9711 

25.8 

7. 0-8.0 

0.9613 

25.8 

8. 0-9.0 

0.9593 

25.8 

9.0-10 

0.9570 

25.8 
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deposited  energy  at  slant  ranges  greater  than  1500  meters.  It  Is  seen 
that  4.3  per  cent  or  less  of  the  source  energy  was  not  deposited  by 
scattering  or  absorption  Interactions  within  1500  meters  slant  range. 
Although  the  maximum  number  of  collisions  allowed  per  history  was  30, 
the  number  actually  undergone  was  less  than  30.  The  average  number 
of  collisions  required  per  history  was  approximately  26  for  source 
energies  above  1.0  MeV  and  reduced  to  13  for  source  energies  between 
0.01  and  0.1  MeV. 

The  neutron  energy  deposited  in  the  atmosphere  out  to  1500  meters 
radius  about  the  source  point  is  shown  in  Table  VI.  Also  shown  is  the  Q 
value  for  neutrons  with  energies  above  2.232  MeV.  The  sum  of  the  energy 
deposited  and  the  Q value  is  also  shown  in  the  table.  The  difference 
between  the  average  source  energy  and  this  sum  represents  the  neutron 
energy  lost  by  neutrons  at  slant  ranges  greater  than  1500  meters  and  by 
the  various  history  termination  methods.  Table  VI  shows  that  no  energy 
was  lost  In  the  neutron  transport  calculations.  The  average  number  of 
collisions  per  history  for  the  neutron  problems  varied  with  the  source 
energy,  with  167  required  for  an  average  source  energy  of  6.93  MeV  to 
107  for  an  average  source  energy  of  0.002175  MeV. 

The  secondary  gamma-ray  energy  generated  per  source  neutron  by 
the  SGSC  procedure  is  given  as  a function  of  the  neutron  source  energy 
Interval  In  Table  VII.  Also  tabulated  Is  the  secondary  gamma-ray  energy 
deposited  within  a slant  range  of  1500  meters.  Only  10  per  cent  of  the 
neutron  collisions  generated  by  the  first  1000  neutron  histories  was 
used  In  generating  secondary  gamma-ray  photons.  An  adjustment  in  the 
weight  of  the  secondary  gamma-ray  source  photons  started  was  used  to 
correct  for  the  bias  sampling.  The  average  number  of  secondary  gamma- 
ray  photon  histories  traced  Is  also  listed  In  Table  VII  for  each  neutron 
source  energy  Interval.  On  the  average,  It  took  approximately  25.6 
collisions  per  secondary  gamma-ray  history  to  deposit  all  of  the  gamma- 
ray  energy.  This  Is  approximately  the  same  number  of  collisions  per 
history  required  for  primary  gamma  rays  with  source  energies  above  2.0 
MeV  (see  Table  V). 

y 


TABLE  VI.  NEUTRON  ENERGY  DEPOSITED  BY  SCATTERING  AND  ABSORPTION,  Q VALUE  AND 
AVERAGE  NUMBER  OF  NEUTRON  COLLISIONS  PER  HISTORY 


Average  Source 
Energy  (MeV) 


Energy  Deposited 
by  Scattering  and 
absorption 


Q Value  Energy  Deposited  Average  No.  of 

Plus  Q Value  collisionji/hl story 


.002175 

.002172 

.00217? 

107 

.006235 

.006281 

.006281 

115 

.01696 

.0168 

.0168 

122 

.0462 

.04607 

.04607 

130 

.1258 

.0258 

.0258 

138 

.2435 

.2425 

.2425 

142 

.4015 

.3997 

.3997 

146 

.6615 

.6610 

.6610 

150 

1.088 

1.089 

1.089 

154 

1.568 

1.546 

1.546 

157 

2.0075 

1.987 

1.987 

159 

2.548'' 

2.538 

.00282 

2.541 

161 

3.2725 

3.250 

.00762 

3.258 

163 

4.8750 

4.694 

.1726 

4.867 

165 

6.930 

5.832 

1.108 

6.940 

167 

8.895 

6.824 

2.063 

8.887 

166 

11.300 

8.390 

2.612 

11.00 

165 

12.75 

9.936 

2.819 

12.76 

165 

14.25 

11.00 

3.270 

14.27 

165 
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Wells  (Ref.  9)  has  published  the  results  of  Monte  Carlo  calcula- 
tions of  the  gamma-ray  dose  rate  in  an  infinite  air  medium  for  point  mono- 
directional  sources  at  source  energies  of  0.6,  1.0,  3.0,  and  6.2  MeV'and 
for  source-receiver  distances  between  500  and  3750  ft.  Zerby  et  al.  (Ref. 
10)  have  also  published  similar  data  for  aourca-raeelvar  distances  of  5, 
10,  20,  40,  65  and  100  ft.  The  gamna-ray  air-scattering  data  from  these 
two  separate  calculations  have  been  Incorporated  as  library  data  in  the 
SKYSHINE  Monte  Carlo  procedure  (Ref.  11)  for  use  in  evaluating  the  gamma- 
ray  exposure  in  air  for  personnel  ip  the  vlclulty  of  the  turbine  building 

of  a nuclear  power  plant.  The  SKYSHINE  library  data  were  transformed  to 

-3  3 

the  air  density  used  in  the  ZAP  calculations  (p-1. 225X10  g/cm  ) for 
comparison  with  the  ZAP  energy  deposition  data.  Figures  3 and  4 show  a 
comparison  of  the  ZAP  calculations  for  6-7  MeV  gamma  rays  at  161  meters 
and  884  meters  with  the  SKYSHINE  data  for  6.2  MeV  gamma  rays.  The  agree- 
ment between  the  calculations  is  good,  considering  the  differences  in  the 
source  energies. 

The  integral  of  the  ZAP-DEPO  calculations  for  line-beam  sources 
over  the  source  polar  angle  6 and  the  azimuthal  angle  (see  Fig.  1)  for 
given  source-receiver  distances  gives  the  energy  deposition  versus  source- 
receiver  distance  for  a point  isotropic  source.  If  ED(E^,0,R)  is  the 
energy  deposited  per  unit  volume  at  position  0,R  for  source  energy  E^, 

Chen 

2lT  IT 

r r , 

ED(E^,R)  " J J ^ ED(E^,0,R)  slnOdOdB 
o 0 

Is  the  energy  daposltad  par  unit  volume  at  position  R for  an  Isotropic 
source  emitting  1 photon  in  all  directions.  Clark  (Raf.  12)  has 
published  the  results  of  moments  method  calculations  of  gamna-ray  snargy 
deposition  in  air.  His  results  wars  reported  in  the  form  of  energy  build- 
up factors  for  a number  of  gaona-ray  energies.  He  expressed  the  build-up 
factor  with  use  of  the  equation 
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D(E^)U(E  )R 

B(vi(E  )R)  = 1 + C(E  )U(E  )R  e ° ° , 

o o o 

where  u(E^)  is  the  total  cross  section  at  energy  E^,  R is  the  source 

receiver  distance  and  E is  the  source  energy.  The  coefficients  C(E  ) 

o o 

and  D(E  ) are  available  from  Ref.  12  for  energies  of  0.5,  1.0,  2.0,  3.0, 
o 

4.0,  6.0,  8.0  and  10  MeV.  These  coefficients  were  plotted  vs  E^  and 
values  were  determined  from  the  plots  at  the  midpoint  energies  of  the 
source  energy  intervals  used  in  the  ZAP-DEPO  calculations.  The  total 
energy  deposition  at  position  R for  a unit  isotropic  source  at  energy 
E^  is  given  by 

-U(E  )R  , 

ED(IJ(E  )R)  = E u (E  )B(u(E  )R)  e ° (AirpR^)  , 
o o ea  o o 

where  macroscopic  energy  absorption  coefficient  and  p 

is  the  air  density.  Values  of  ED(u(E^)R)  were  computed  using  values 

of  u (e  )/p  from  Ref.  13  and  R values  given  by 
ea  o 

for  the  radial  Intervals  R^  to  R^  listed  in  Table  II  and  and  source 
energies  E^  of  4.5,  8.5,  7.5,  6.5,  5.5,  4.5,  3.5,  2.5,  1.5,  and  0.75. 

Table  VIII  lists  the  ratio  of  the  ZAP-DEPO  calculations  of  the  energy 
deposition  to  the  Clark  results.  Also  shown  in  Table  VIII  are  the  results 
of  integrating  the  energy  deposition  over  the  volume  extending  out  to 
1500  meters.  The  poorest  agreement  between  the  ZAP-DEPO  data  and  Clark's 
data  is  at  the  first  slant  range  for  all  source  energies  and  at  the  last 
slant  range  for  source  energies  less  than  3.5  MeV.  In  general,  the 
agreement  is  good,  less  than  10%  difference,  for  most  slant  ranges  and 
source  energies.  The  ZAP-DEPO  data  tends  to  slightly  overpredict  the 
Clark  data.  The  ratios  of  the  total  energy  deposition  in  a volume  out 
to  a radius  of  1500  meters  shows  that  the  ZAP-DEPO  data  give  about  5 
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TABLE  VIII.  RATIO  OF  THE  GAMMA-RAY  ENERGY  DEPOSITION  PER  UNIT  VOLUME  OF  AIR 
FROM  A POINT  ISOTROPIC  SOURCE  AS  COMPUTED  WITH  ZAP-DEPO  TO  THAI 
COMPUTED  USING  CLARK'S  ENERGY  BUILDUP  FACTOR  DATA  VERSUS  SOURCE 
ENERGY  AND  SLANT  RANGE 


SOURCE  ENERGY  (MeV) 


R(cm)  1 

EBa 

1.5 

2.5 

3.5 

4.5 

5.5 

6.5 

7.5 

8.5 

7.937+2* 

1.885 

1.967 

1.823 

1.929 

1.817 

1.949 

1.889 

1.902 

1.831 

7.0^0 

1.651+3 

1.081 

1.154 

1.151 

1.111 

1.033 

1 . 230 

1.167 

1.123 

1.037 

1.260 

3.302+3 

1.088 

1.165 

1.146 

1.173 

1.160 

1.185 

1.195 

1.171 

1.180 

’ .211 

5.192+3 

1.035 

1.048 

1.072 

1.034 

1.049 

1.063 

1.060 

1.055 

1.143 

1 . Ill 

7.140+3 

0.938 

1.035 

1.112 

1.055 

1.038 

1.103 

1.147 

1.004 

0.984 

1.140 

9.110+3 

0.946 

0.939 

0.997 

1.024 

1.061 

1.006 

1.038 

1.017 

1.010 

1.031 

1 . 139+4 

0.976 

0.972 

0.971 

0.949 

0.974 

1.056 

1.070 

1.085 

1.089 

i.o«; 

1 . 386+4 

0.985 

0.977 

0.995 

1 . 064 

1.075 

1.072 

1.025 

0.971 

1.117 

1.028 

1.635+4 

0.996 

0.985 

0.955 

1.092 

1.042 

1.036 

0.981 

1.088 

1.068 

1.075 

1.883+4 

0.987 

0.993 

1.020 

1.106 

1.106 

1.054 

0.980 

1.030 

1.069 

1.049 

2.277+4 

1.002 

1.040 

1.051 

1.056 

1.048 

0.986 

1.071 

1.088 

1.071 

1.01' 

2.773+4 

1.061 

1.021 

1.012 

1.001 

0.780 

1.089 

1.024 

1.073 

1.056 

I.0S2 

3.269+4 

1.142 

1.095 

1.067 

1.053 

1.042 

0.994 

1.010 

1.072 

1.081 

1 .042 

3.767+4 

1.133 

1.051 

1.030 

1.005 

1.020 

1.018 

0.953 

1.000 

1.058 

1.033 

4.265+4 

1.205 

1.034 

1.079 

1.057 

1.077 

1.019 

1.017 

1.039 

0.989 

0.9  '9 

4.763+4 

1.286 

1.043 

1.057 

1.076 

1.114 

1.142 

1.088 

1.016 

1.022 

1 . 04  1 

5.545+4 

1.354 

1.092 

1.118 

1.112 

1.079 

1.062 

1.028 

1.045 

1.004 

1.027 

7.140+4 

1.548 

1.271 

1.279 

1.143 

1.108 

1.075 

1.115 

1.071 

1.011 

1.04  4 

9.110+4 

1.375 

1.208 

1.163 

1.086 

1.080 

1.005 

1.059 

1.044 

0.993 

0.968 

1.109+5 

1.166 

1.199 

1.125 

1.077 

1.048 

1.047 

1.114 

0.932 

1.026 

O.Ol'i 

1.366+5 

2.487 

1.588 

1.477 

1.174 

1.039 

1.130 

1.028 

1.063 

1.095 

0.913 

0-1. 5+5 

1.070 

1.066 

1.082 

1.078 

1.073 

1.068 

1.066 

1.092 

1.056 

1.053 

*Read  7.937+2  as  7.937X10^ 
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to  8%  more  energy  deposition  than  do  the  Clark  data. 

There  does  not  appear  to  have  been  any  previous  calculations  of 
the  neutron  and  secondary  gamma-ray  dose  rates  in  air  for  line-beam 
neutron  sources  other  than  those  reported  previously  by  Wells  in  Ref.  14. 
Those  calculations  were  performed  in  1960  and  the  neutron  and  secondary 
gamma-ray  production  cross  sections  used  in  those  calculations  are  out  of 
date.  In  addition,  the  calculations  were  not  made  for  source-receiver 
distances  greater  than  100  ft.  Therefore,  the  ZAPN-DEPO  and  ZAPGAM-DEPO 
calculations  for  line-beam  sources  were  integrated  over  the  source  polar 
angle  to  provide  data  for  comparison  with  calculations  for  point-isotropic 
neutron  sources.  Straker  (Ref.  15)  has  reported  the  results  of  his  cal- 
culations of  the  neutron  dose  (air  kerma)  in  an  infinite  homogeneous  air 

medium  for  a 12.2-15  MeV  point-isotropic  neutron  source.  Straker's  data, 

2 2 

in  units  of  4tiR  times  the  neutron  dose  (cm  Rad/source  neutron),  are 

compared  in  Fig.  5 with  the  results  given  by  the  ZAPN-DEPO  calculations 

for  a neutron  source  emitting  uniformly  between  12.0  MeV  and  15  MeV. 

-3 

The  dose  rates  computed  by  ZAPN-DEPO  for  an  air  density  of  1.225X10 

3 -3  3 

g/cra  were  transformed  to  an  air  density  of  1.11X10  g/cm  , which  was 

the  air  density  used  in  Straker's  calculations.  The  comparison  in  Fig.  5 
shows  that  the  ZAPN-DEPO  calculations  slightly  overpredict  Straker’s  cal- 
culations as  a function  of  slant  range  by  about  10  15%  over  most  of 

the  slant  range  considered.  Considering  the  possible  differences  in  the 
neutron  cross  sections  used  in  the  two  different  calculations,  the  agree- 
ment between  the  calculations  is  good. 

Straker  (Ref.  15)  has  also  reported  the  results  of  calculations 
of  the  secondary  gamma-ray  dose  in  an  infinite  homogeneous  air  medium 
that  results  from  the  atmospheric  absorption  and  scattering  of  neutrons 
that  were  emitted  by  a 12.2  MeV  to  15  MeV  point-isotropic  source.  His 
results  are  compared  in  Fig.  6 with  similar  results  from  the  ZAPGAM- 
DEPO  calculations  for  a 12.2  to  15  MeV  point  isotropic  neutron  source. 
Straker's  calculations  were  based  on  the  use  of  Young  and  Foster’s 
(Ref.  16)  evaluation  of  the  neutron  cross  sections  and  the  secondary  gamma- 


RANGE  (l 


f ZAPN  Calculat 
'a  (Ref.  15)  Ca 


Fig.  6.  Conparlson  of  ZAPGAM  Calculations  of  Secondary  Gamma-Ray  Dose  In  Air 
with  Data  from  Straker  (Ref.  IS)  Calculations 
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ray  production  cross  sections  for  nitrogen.  Young  and  Foster's  eval- 
uation was  to  later  become  the  basis  of  the  ENDF/B-IV  evaluated  neutron 
cross  sectlon.s  for  nitrogen.  Straker's  data  are  about  35X  higher  than 
the  ZAPGAM-DEPO  data  for  slant  ranges  less  than  1500  meters.  The 
differences  between  the  two  calculations  are  believed  to  be  due  to 
differences  in  the  secondary  gamma-ray  production  cross  sections  'jstj 
in  the  two  calculations. 

The  DEPO  calculations  sum  the  energy  deposltlcn  at  each  cullisivii , 
with  the  time  dependence  given  in  terms  of  retarded  time  where  retarded 
time  is  the  actual  time  at  collision  minus  the  time  it  would  take  for 
light  to  propagate  from  the  source  to  the  collision. 

The  time  dependence  of  the  gamma-ray  energy  deposited  in  a 

sphere  of  air  1500  meters  in  radius  per  keV  of  source  energy  for  point- 

isotroplc  sources  emitting  within  each  of  the  12  source  energy  intervals 

is  given  in  Table  IX.  The  gamma-ray  energy  deposited  in  the  atmosp^-er? 

-20 

within  a retarded  time  interval  of  0 to  10  seconds  represents  the 

energy  deposited  by  first-order  scattering  and  absorption  interactluna . 

It  is  seen  that  the  fraction  of  the  source  energy  that  is  absorbed  as  a 

result  of  first-order  scattering  or  absorption  events  reduces  with 

decreasing  source  energy  until  a source  energy  of  10  keV  is  reached  and 

then  the  fraction  Increases  with  decreasing  source  energy.  The  fraction 

of  the  gamma-ray  source  energy  absorbed  within  a sphere  of  air  with  a 

-1  -20 

1500-meter  radius  (keV  sec  /keV  of  source  energy  times  10  “ sec)  closely 

follows  the  variation  of  the  ratio  of  the  energy  absorption  to  total 

cross  section  with  energy  (see  Fig.  7),  where  the  cross  section  data  were 

taken  from  Ref.  13.  The  time  dependence  of  the  absorbed  gamma  radiation 
-20  -5 

for  times  between  10  sec  and  1X10  sec  reduces  rapidly  with  an  increase 
in  retarded  time.  The  rate  of  reduction  with  increasing  retarded  t Ime 
reduces  with  a decrease  in  the  source  energy.  At  a retarded  time  of 
7X10  ^ sec  there  is  a factor  of  10  more  energy  being  deposited  by  the 
lO-to-lOO-keV  source  than  is  being  deposited  by  the  9-to-lO-keV  source. 
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The  time  and  radial  distance  dependence  of  the  energy  deposition 
-3  -1 

(keV  m sec  /keV  of  source  energy)  from  a 6-to-7-MeV  point-isotropic 

gamma-ray  source  is  shovm  in  Table  X.  The  energy  deposition  rate  in 

“8 

retarded  time  reduces  with  increasing  time  for  times  after  2X10 

seconds  at  about  the  same  rate  for  each  radial  Interval.  Figure  8 shows 

a plot  of  the  energy  deposition  rate  for  radial  intervals  of  175  to 

200  meters,  500  to  600  meters  and  1000  to  1200  meters  for  retarded  times 
—8  —5 

between  10  sec  and  10  sec.  The  energy  deposition  rate  varies 

“1  33  ”8 

approximately  as  a power  curve  with  ED(t,R)»a(R)t  ' for  t>2X10  sec 
and  radial  distances  R^1500  meters.  The  coefficient  a(R)  varies  with  the 
radial  distance. 

The  time  distribution  of  the  neutron  energy  deposition  in  a sphere 
of  air  1500  meters  in  radius  is  shown  in  Table  XI  as  a function  of  the 
source-energy  Interval.  In  the  earlier  time  Intervals  the  energy  deposition 
rate  Increases  by  a factor  of  approximately  2X10^  as  the  source  energy 
interval  Increases  from  0.001  to  0.00335  MeV  to  13.5  to  15  MeV.  With  an 
increase  in  retarded  time  the  energy  deposition  rate  Increases  slower 
with  an  Increase  in  the  source  energy.  At  retarded  times  after  2X10  ^ 
sec,  the  energy  deposition  rate  Increases  with  source  energy  to  a maximum 
after  which  it  decreases  with  a further  increase  in  the  source  energy. 

The  source  energy  at  which  the  maximum  energy  deposition  rate  is  found 
decreases  with  increasing  retarded  time.  The  rate  at  which  the  energy 
deposition  reduces  with  increasing  retarded  time  increases  with  increasing 
source  energy. 

The  secondary  gamma-ray  energy  deposition  as  a function  of  tisw* 
and  neutron  source  energy  interval  in  a 1500-meter  sphere  about  a point- 
isotropic  neutron  source  is  given  in  Table  XII.  For  source  energies  less 
than  3.68  MeV  all  of  the  secondary  gamma  rays  are  produced  by  n,Y  inter- 
actions. Since  the  n,Y  cross  section  Increases  with  reducing  neutron 
energy,  all  of  the  capture  gamma  rays  are  produced  at  late  times.  A 
significant  portion  of  the  total  secondary  gamma  emission  occurs  in  the 


I 


TABLE  X.  TIME-DEPENDENT  ENERGY  DEPOSITION  IN  AIR  VS 

RADIAL  DISTANCE  FROM  A POINT  ISOTROPIC  GAMMA- 
RAY  SOURCE  EMITTING  UNIFORMLY  IN  THE  ENERGY 
INTERVAL  FROM  6 to  7 MeV 
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100  to  1200  METER 
RADIAL  DIST/uNCE 


RETARDED  TIME  (sec) 


Time  Dependence  of  the  Energy  Deposition  from  a 6 to  7 
MeV  Point  Isotropic  Gamma-Ray  Source 
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-2  -1 

last  time  interval  (1.5X10  to  1.0X10  sec)  when  the  neutron  source 
energy  is  less  than  3.68  MeV,  When  the  source  energy  is  above  3.68  MeV, 
most  of  the  secondary  gamma-ray  energy  deposition  comes  from  inelastic 
gammas  and  is  produced  by  neutrons  that  have  undergone  Inelastic  colJi- 
sions  before  about  5X10  ^ seconds.  The  total  secondary  gamma-ray  energy 
deposited  within  the  1500  m radius  decreases  slowly  with  an  increase  in 
the  neutron  source  energy  until  the  neutron  source  energy  exceeds  3.68 
MeV  after  which  it  increases  to  a maximum  for  neutrons  in  the  12. G to 
13.5  MeV  interval  and  then  reduces  with  increasing  energy. 

Of  interest  to  a number  of  organizations  working  on  problems 
relating  to  exposure  of  personnel  and  equipment  to  the  radiations  pro- 
duced by  nuclear  sources  is  the  time-dependent  energy  deposition  vs 
range  in  air  for  point  isotropic  sources  and  the  time-independent  energy 
deposition  vs  range  for  line-beam  sources.  The  time-dependent  energy 
depositions  in  the  21  radial  Intervals  between  0 and  1500  meters  listed 
in  Table  II  for  point  Isotropic  nuclear  sources  emitting  in  the  discrete- 
energy  Intervals  listed  in  Table  IV  are  tabulated  in  Vols.  II,  III,  and 
IV  for  primary  gammas,  neutrons  and  neutron-induced  secondary  gamm.i  rays 
respectively.  The  energy  deposition  data  given  in  Vols.  II  through  IV 
have  the  following  units: 


Vol.  No. 

Source 

Units  of 

Energy  Deposition  Data 

II 

gamma  rays 

keV  m ^ 
-3 

sec  ^/keV  of  source  energy 
-1, 

III 

neutrons 

keV  m 

-3 

sec  /source  neutron 
-1, 

IV 

secondary 
gamma  rays 

keV  m 

sec  /source  neutron 

The  radial  distances  are  in  meters  and  the  retarded  times  are  In  seconds. 

-20 

For  the  gamma-ray  sources,  the  energy  deposition  in  the  O-to-10  -second 
delay-time  interval  is  taken  to  be  produced  by  the  first-order  scattering 
and  absorption  interactions  undergone  by  the  direct  radiation.  The  gamnui- 
ray  energy  deposition  computed  for  retarded  times  greater  than  10  seconds 
in  a given  radial  interval  is  that  produced  by  gamma  rays  that  underwent 
their  first-order  collision  outside  of  the  deposition  volume  and  then  iinder- 
went  second-order  or  greater  collision  within  the  deposition  volume. 
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The  time-integrated  total  energy  deposition  data  vs  radial 
distance  and  source-energy  Interval  are  also  listed  in  Vols.  II, 

III  and  IV  for  point-isotropic  gamma-ray,  neutron  and  neutron- 
produced  secondary  gamma-ray  sources,  respectively.  The  units  for 

_3 

the  gamma-ray  energy  deposition  data  are  keV  m /keV  of  source  energy. 
The  units  for  the  neutron  and  secondary  gamma-ray  energy  deposition 

_3 

data  are  keV  m /source  neutron. 

The  time- Integrated  energy  deposition  data  for  line-beam  sources 
as  a function  of  radial  distance,  deposition  polar-angle  interval 
(Table  II)  and  source-energy  interval  have  been  curve  fitted  by  a 
least  squares  method  with  the  use  of  the  equation 


-A,  , B -C^  , R, 

ED(9,,R.,E,  ) - e e ^ 

1 j k j 

where  1 and  k are  indices  denoting  the  deposition  polar  angle  interval 
and  source  energy  Interval  and 


R 


i 


+ R 


J 


for  R and  R being  the  lower  and  upper  bounds  of  the  jth  radial  interval. 

3“-^  J 2 

The  parameters  A,  B and  C and  the  coefficient  of  determination  r as  a 

function  of  the  deposition  polar  angle  and  source  energy  Intervals  are 

2 

given  in  Vol.  V.  Values  of  r close  to  1.00  Indicate  the  goodness  of  the 

fit  to  the  DEPO-calculated  energy  deposition  data.  The  energy  deposition, 

— —3 

ED(9  ,R  ,E  ) is  in  units  of  keV  m /keV  of  source  energy  for  gamma-ray 

i J _3 

sources  and  keV  m /source  neutron  for  neutron  and  secondary  gamma-ray 
sources . 


The  time-integrated  energy  deposition  data  given  in  Vols.  II  through 

-3 

V can  be  converted  to  air  kerma  rate  (exposure  rate  in  air,  p “ 1.225X10 
g/cm^)  by  multiplying  the  energy  deposition  data  with  one  of  the  following 


factors: 
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Source 

Kerma  Factor 

Gamma  Rays 

4.708X10”^^ 

Neutrons 

4.708X10'^^ 

Secondary  Gamma  Rays 

4.708X10”^^ 

The  units  of  the  resulting  air-kerraa  rate  is  rad  hr"^/source  particle 
sec  and  is  the  average  source  energy  within  the  source  energy 
interval. 
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III.  ENERGY  DEPOSITION  DATA  FOR  CONICAL  SOURCES 


A modified  version  of  the  DEPO  program  was  developed  for  use  in 
determining  the  spatial  and  time  distributions  of  the  nuclear  energy 
deposited  in  an  infinite  air  medium  for  conical  sources  from  the  data 
generated  for  line-beam  sources.  For  each  point  of  energy  deposition 
recorded  on  a ZAP,  ZAPN  or  ZAPGAM  history  tape  for  each  of  the  mono- 
directional  sources,  a cosine  value  between  the  limits  of  the  conical 
source  bounds  was  selected  at  random  (see  1**  Fig*  !)•  The  new  direction 
selected  is  considered  to  be  the  new  source  direction  and  the  polar  dis- 
placement angle,  0 , from  the  monodirectional  source  was  referenced  to  this 
E 

new  source  direction.  A number  of  random  azimuthal  angles,  y,  about  the 
new  source  direction  was  selected  uniformly  between  0 and  Ztt.  The  number 
of  azimuthal  angles  selected  was  established  as 

N - 2*IFlX(^rccose„)  + 3 

where  a was  slightly  less  than  5.  The  number  of  points  of  energy  deposition 
(each  depositing  ^ of  the  original  energy  deposition)  ranged  between  3 and 
11.  The  third  side  of  the  spherical  triangle,  6,  references  the  points  of 
energy  deposition  back  to  the  axis  of  the  conical  source.  The  angular 
position  of  the  energy  deposition  for  a conical  source  was  recorded  in 
te^ms  of  the  angle  6,  the  angle  between  the  point  of  deposition  and  the 
axis  of  the  conical  source.  The  radial  distance  and  the  time  to  collision 
was  taken  to  be  the  same  as  that  computed  for  the  monodirectional  source. 

The  bounds  of  the  source  polar  angle  Intervals  for  the  conical 

sources  are  0-5,  5-10,  10-20,  20-30,  30-40,  40-50,  50-60,  60-70,  70-80, 

and  80-90  degrees.  Deposition  data  were  generated  with  the  revised 

version  of  DEPO  for  conical  source  polar  angles  up  to  90°.  Energy 

deposition  data  for  conical  source  polar  angles  greater  than  90°  can 

be  obtained  by  reordering  the  polar  angle  intervals  in  the  deposition 

geometry.  The  new  deposition  polar  angle  bounds  0^  for  k ■ 0,  21  are 

defined  by  6'  » 180°  - 0.,  , when  90  < 0 < 180°. 

k 21-k  s — 
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Energy  deposition  data  for  source  polar-angle  bounds,  4/^,  for 
a conical  source  between  100*^  and  110*’  are  the  same  data  as  those 
computed  for  conical  source  polar-angle  bounds  of  70^  and  80°  wl th  the 
deposition  polar  angles  0^  for  the  100°  to  110°  conical  source  defined 


as  given  above  In  terms  of  the  deposition  polar  angles  0^^  for  the 
70°  to  80°  conical  sources. 


i 


1 
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IV.  RENDER  PROCEDURE 

A computer  procedure,  designated  as  RENDER,  was  developed  for  use 
in  computing  the  time  and  spatial  dependent  nuclear  energy  deposition  In 
the  atmosphere  for  time,  energy  and  direction  dependent  point  sources. 

The  energy  deposition  at  spatial  position  9^^  ^ 6 ^ ®i+l 
tj  ^ r £ ’^j+l  time  T Is  given  by  the  Integral  equation 


1 01+1.  rj  1 r < rj^^.  T) 


I r r-'  r 

E (ti“0  0^^  r^  il^-O  0 


j+1 


S(E,iJ;,t)  T(E,0,r,tJ;,T-t)slniJ)r  sln0dtdi^drd0d(j)dE 


(1) 


I®  1+1  l’^j+1 

■0  0^  tj  \|>“0 


sinter  sln0dtWrd0d(|) 


where  S(E,ip,t)  Is  the  energy,  direction  and  time-dependent  source  term 

T(E,0,r,i/),T-t)  Is  the  energy  deposited  per  unit  volume  for  Initial 

source  energy  E emitted  In  direction  ijj  and  deposition 
occurring  at  position  r,0  and  time  T-t. 

The  convolution  Integral  portion  of  Equation  (1) 

I S(E,iiJ,t)  T(E,0,r,il),T-t)dt  (2) 

0 

Is  evaluated  in  RENDER  by  assuming  analytical  functions  for  the  source 
S(E,i(),c)  and  the  transmission  T(T-t).  Ignoring  the  dependence  on  source 
conical  angle  ip,  source  energy  E and  spatial  location  r,0  we  can  rewrite 
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the  above  equation  as 


R(T)  - 


S(t)  T(T-t)dt 


(3) 


where  S(t)  and  T(T-t)  are  analytic  functions.  In  practice,  however, 
source  values  are  given  at  discrete  time  points,  say  transmission 
data  Is  available  as  average  arrival  rates  between  two  retarded  times, 
say  X j ■ Considering  these  partical  constraints  and  the  fact 
that  the  convolution  resultant  R Is  to  be  defined  at  retarded  times  r, 
a transformation  of  variables  was  made  In  Eq.  3 as 

X • T-t  , 


thus 


R(t)  - J T(x)  S(T-x)dx  . 
0 


(4) 


Since  T(x)  Is  given  as  a constant  (average)  between  time  limits,  say 
T (Xj),  Eq,  4 can  be  written  as 


R(T)  - T(x^) 


S(T-x)dx  + • • • + T(x 


.>r 


S(T-x)dx 


J-1 


U 

-Jtcxj) 


j-1 


mln(T,x,) 


S(T-x)dx 


J-1 


*0  ~ ® Index  of  the  Inequality 

*J-1  ^ • 

Reversing  the  previous  change  of  variables  gives 


(5) 


U « T-X 


R(T)  - y T(Xj] 

J-1 


J-1 

S(u)du 

T-min(r,Xj) 


The  integral  part  of  Eq.  6 could  be  solved  in  a number  of  ways.  Although 
the  form  of  S was  assumed  to  be 

S(u)  - aue  , (7) 

the  RENDER  program  was  written  such  that  only  minor  programming  changes 
would  be  necessary  if  another  form  of  S were  considered  to  be  more 
applicable.  Coefficients  a and  b in  Eq.  5 are  regenerated  for  each 
pair  of  source  times,  u^  and  as 

b - In  (— 1 

“i+l~"i  Ti+l  ^i  I 


i b u. 
a * — e 1 


These  coefficients  cause  S(u)  to  be 


- a u^  e ** 


S,  . , ” a u 


1+1  i+1 


Results  of  the  Integral  portion  of  Eq.  6 follow  as 


I S(u)du  - ^ 


e“*’“  (bu+1)  . 
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Convolution  results  can  thus  be  obtained  as 

T I, 


R(t) 


J-i 


i-i. 


■ -bU*  -bU 

e (bUi^^+1)  - e (bU2  ^1-1) 


(8) 


where  Ij^  and  I2  are  developed  from  the  following  Inequalities 


u < T - min  (t,  X ) < u 
1 ^ ^1 


Ut  < t - X,  , < u 


J~1  - -I,+l 


and 


u =0 
o 


Values  of  U|  ^ and  U2  ^ are  give  by 


- max  (T-mln  (t.x^), 

'^2,1"  ^'^“*J-1’  '*12+2-1^ 

and  coefficients  a and  b are  generated  at  each  new  value  of  1,  with 
the  exception  of  1 ■ 0,  The  form  of  the  source  term  In  Eq.  7,  considered 
on  the  Interval  0 £ u < “>,  has  some  Important  properties.  First 


S(u)  - 0 when  u - ^ 


d , 1.  ab 

du^  ** 


which  for  a>0,  b>0  gives  a single  maximum  et  u-^ 
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Secondly 


11m  c(  \ au 

u-*“>  U"*®°  bu 

e 


lim  .0.0 

U-Kd  d bu»  ® 


and 


S(u“0)  “ a"0*e  ^^0)  _ q 


Thus  the  form  of  the  source  has  the  correct  boundary  conditions  and 
contains  one  maximum. 

To  test  the  effectiveness  of  the  numerical  solution  of  the 
convolution  process,  an  analytic  source  and  transmission  were  chosen 

, ,-12  -lO^u 

S(u)  ‘10  u e 


T(x)  - 2X10^  * 


which  gives 

R(T)  - 3.125X10^  ^ {l  - ^ I8Xlo\+l]}  . j 

Graphs  of  S,  T and  R are  shown  In  Fig.  9.  Since  the  form  of  the  source 
was  Identical  to  the  assumed  form  of  the  source.  It  would  seem  reasonable 
that  the  results  of  the  numerical  solution  be  Independent  of  the  range  of 
the  source  or  the  number  of  discrete  values  used  (minimum  of  two  are  neces- 
sary to  evaluate  the  coefficients  a and  b) . This  was  Indeed  the  case. 

The  only  errors  possible  were  In  1)  computer  accuracy  and  2)  limiting  values 
used  to  generate  the  average  transmission  value  from 
rx. 


i 


T(x)dx 


T(x^)  - 


-2X10"x 


J-1 


-2X10  X 


j 


“J 

dx 

J-1 


*J  ■ *j-l 


iawiiwiiiimiiBH'i 


iitimi! 


taiiiliHIHU!;; 

IlliHDpffili 


inirmi: 


3.125X10 


TIME  (sec) 

Analytic  Functions  used  In  Test  Problem 


IH 

m 

i| 

sj 

A8 


where  x =0.  Assuming  small  computer  accuracy  errors,  the  range  of 
o 

{x^}  was  varied  but  the  number  of  x^  was  held  constant  at  30.  Figure  10 
gives  a tabulation  of  the  per  cent  error  in  the  RENDER  results  of  five 
numerical  solutions  of  the  test  problem.  In  each  case  the  transmission 
time  bounds  were  calculated  as 


X 


j 


XeJ*" 


where  A was  defined  as 


5*10 


X 


X^  was  input  as  10  10  10  10  ^ and  10  This  increase  in  the 

range  of  x^  gave  two  effects;  as  X^  became  smaller  a larger  portion  of 
the  transmission  function  was  defined  but  each  value  of  the  transmission 
became  less  accurate  due  to  the  Increased  time  coverage.  Although  the 
transmission  time  set  {x^}  will,  in  practice,  be  fixed,  it  la  evident 
that  the  X^  should  be  small  enough  to  Include  all  contributions  to  the 
transmission  at  early  times  and  as  many  values  of  x^  as  possible  within 
computer  mamory  storage  requirements. 

A second  sample  RENDER  problem  was  run  using  a source  term  and 
transmission  function  that  was  dependent  on  the  source  polar  angle  ijj. 

The  source  term  was  defined  by 


S((^,t)  " (l+cosii')te  (source  units  per  steradian) 

ar.d  the  transmission  function  was  defined  by 

T(0,\J(,r,t,E)  - — j — I e ^^'slnGslniJ;  . (9) 

r (r>l) 

A value  of  y “ A55. 8866187  was  used  and  symmetry  in  angle  for  T (le. 
T(it-0 ,17-11;)  - T(0,4/)  was  employed  to  ease  the  analytic  solution.  The 
failure  of  the  sample  problem  to  typify  the  type  of  production  problems 
expected  for  RENDER  lies  mainly  in  the  two  angle  terms : 0 and  ij;  . A 
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Percent  Error 


RETARDED  TIME 

Lower 

Bound  of 

Xo(sec) 

T(sec) 

lO"^ 

1 

o 

1 

o 

lO"® 

•- 

1 

1.  lOOOE-07 

-58.04 

-9.91 

-.46 

-.01 

-.01 

1 . 2000E-07 

-50.02 

-9.8  \ 

-.41 

-.01 

- .02 

1. 4000E-07 

-57.96 

-9.72 

-.30 

- . 01 

-.  02 

1 . 6000E-07 

-57.90 

-9. 59 

-.22 

-.02 

-.02 

1.0OOOF-O7 

-57.04 

-9.47 

-.18 

- .02 

- .03 

2.1000E-07 

-57.75 

-9.27 

-.13 

-.02 

-.03 

2.3000E-07 

-57 . 69 

-9. 14 

-.11 

- . 02 

-.03 

2. 7000E-07 

-57.57 

-8.88 

-.00 

-.03 

- .04 

3.  OOOOE-07 

-57. 48 

-8.68 

-.06 

- .03 

- . 05 

3.4000E-07 

-57.35 

-0.42 

-.05 

-.03 

-.04 

3.9000E-07 

-57. 19 

-8.08 

-.05 

-.04 

-.06 

4 .4000E-07 

-57.03 

-7.73 

-.04 

-.04 

-.07 

5. OOOOE-07 

-56.83 

-7.31 

- . 04 

-.05 

- . 07 

5.600nE-07 

-56 . 63 

-6.07 

-.04 

-.05 

- .09 

6 . 400nE-07 

-56.36 

-6.28 

-.04 

-.07 

- . 10 

7 . 4000E-07 

-56.00 

-5.51 

-.04 

-.06 

-.10 

8 . 2000E-07 

-55.71 

-4.80 

-.04 

- .08 

-.14 

9.4000E-07 

-55.25 

-3.90 

-.05 

-.09 

-.14 

1.000f)E-06 

-55.01 

-3.40 

-.05 

-.06 

- .13 

1 .2000E-06 

-54 . 19 

-1.69 

-.06 

- . 13 

-.22 

1 . 3000E-06 

-53.76 

-1 . 08 

-.06 

-.11 

-.2? 

1 .5000E-06 

-52.85 

-.32 

-.07 

-.17 

- .21 

1 . 7000E-06 

-51.87 

. 09 

-.10 

- . 14 

-.33 

1 .9000E-06 

-50.83 

. 29 

- . 09 

-.20 

-.32 

2. 2000E-06 

-49.14 

. 43 

-.  13 

-.21 

33 

2.5000E-06 

-47.20 

. 45 

-.15 

-.27 

-.53 

2. 8000E-06 

-45.23 

. 42 

-.  13 

-.34 

- . 38 

3.2000E-06 

-42.20 

.35 

-.22 

-.29 

- . 53 

3.6000E-06 

-30.80 

.27 

-.17 

-.50 

-.05 

4.20D0E-06 

-32.92 

. 18 

-.18 

-.31 

-.30 

4.0OOnE*O6 

-26.02 

.10 

-.36 

- .72 

-1.  10 

5.4000E-06 

-17.98 

.01 

-.17 

-.  16 

-1.20 

6. 0000E*06 

-8. 69 

-.07 

-.46 

- . 09 

. 1 1 

6.8000E-06 

5.90 

-.  14 

-.26 

-.  32 

-1.66 

7.8000E-06 

28.19 

-.19 

.01 

-1.00 

-1.39 

0. OOOOE-06 

55.78 

-.24 

-.77 

-.  75 

1.31 

1.  OOOOE'05 

97.43 

-.27 

.44 

- .19 

-3.55 

1 . lOOOE-05 

127.72 

-.00 

-1  . 00 

-2.35 

-1.87 

1 . 200CE-05 

101.  19 

..^2 

1.51 

1.97 

6.64 

1 .4000E-05 

38 .16 

. 61 

-.32 

-3. 14 

-5 .45 

1.6000E-05 

11.25 

.89 

-1.32 

7.13 

2. 92 

1.8000E-05 

3.0  7 

1 . 75 

6 . 06 

-2.51 

13.37 

2. OOOOE-05 

. 03 

1.45 

-2.  73 

.01 

-lO.  Dl 

2.3000E-05 

.22 

3.66 

5.77 

5.55 

16.67 

2. 6000E-05 

.JO 

3.99 

4.35 

-1.50 

20.86 

3 . OOOOE-05 

.36 

0. 45 

-5.65 

19.22 

-0 . 35 

3. 4000E-05 

. 21 

11.63 

26.96 

-1.00 

85.  76 

3.0OOOE>O3 

.56 

7.01 

-.65 

92.22 

4.87 

4. 4000E-05 

. 21 

22.55 

-7.06 

-8 .03 

26.47 

5.0000E-05 

2.04 

29.67 

80.67 

100.04 

315 . 42 

Fig.  10.  Percent  Error  In  RENDER  Calculation  for  Sample  Problem 
As  a Function  of  the  Lower  Bound  of  X 
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transmission  term  which  Is  separable  In  6 and  \p  Is  obviously  not  repre- 
sentative of  nature  and  a simple  sine  dependence  does  not  reflect  the 
variations  found  In  the  transmission  data  which  has  been  previously 
generated.  A more  suitable  dependence  on  the  source  and  deposition 
angle  would  possibly  have  been  a function  of  the  form 

1(6,4^)  = exp  [k|sln'^((|;-9)  I ] 

Use  of  such  a function  would,  however,  complicate  the  analytic  solution 
to  such  a degree  that  numerical  methods  would  be  required  to  obtain  results. 
This  would  negate  the  purpose  of  the  sample  problem.  Employment  of  the 
simple  function  (Eq,  9)  will  produce  errors,  but  predictable  errors  which 
can  be  used  to  both  demonstrate  Inherent  errors  of  the  procedure  as  well 
as  the  accuracy  of  the  procedure  to  generate  the  .»pproprlate  solution. 

The  desired  solution  Is 

R(0i  1 0 < 0^^^.  r.  < r < rj^j,  X,  E^  < E < E^^^  - 

3tt  ^^l+r^l  - ccse^^^,sln9.^^)(tan~^^^^  - tan  ^ 

^ (cosO.  - cos0.^^)(r]^^  - rj^) 


{e'^"[(Y-a)T-Hl  H-  e~^") 

(Y-a)^ 

To  allow  the  source  and  transmission  to  be  Input  to  RENDER,  bounds  were 

established  for  all  variables.  There  were  34  equally  spaced  source  angle 

Intervals  {ti  such  that  0 < i{/  < tt  and  20  source  times  t equally  spaced  In 

—8  —5 

logarithm  such  that  10  £ t £ 10  seconds.  Transmission  data  was  defined 

for  17  source  angles  such  that  0 ^ 5.  7 » times  t such  that 

-8  -2  ^ 

10  < t < 10  , radial  and  deposition  angles  as  shown  In  Table  XIII. 

Both  source  and  transmission  were  Independent  of  energy. 
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TABLE  XIII.  SAMPLE  PROBLEM  PARAMETERS 


nterval 

Number 

Output 

Polar  Angle 
Interval 
(degrees) 

Radial 

Interval 

(meters) 

Deposition 
Polar  Angle 
Interval 
(degrees) 

1 

0 

1 

0 - 

10 

0 

1 

2 

1 

- 

2 

10  - 

20 

1 

- 

2 

3 

2 

- 

5 

20  - 

40 

2 

- 

3 

4 

5 

- 

10 

40  - 

60 

3 

- 

5 

5 

10 

- 

20 

60  - 

80 

5 

- 

7 

6 

20 

- 

30 

80  - 

100 

7 

- 

10 

7 

30 

- 

40 

100  - 

125 

10 

- 

15 

8 

40 

- 

60 

125  - 

150 

15 

- 

20 

9 

60 

- 

70 

150  - 

175 

20 

- 

30 

10 

70 

- 

80 

175  - 

200 

30 

- 

40 

11 

80 

- 

90 

40 

- 

50 

12 

90 

- 

100 

50 

- 

60 

13 

100 

- 

110 

60 

- 

70 

14 

110 

- 

120 

70 

- 

80 

15 

120 

- 

130 

80 

- 

90 

16 

130 

- 

150 

90 

- 

100 

17 

150 

- 

160 

100 

- 

120 

18 

160 

- 

170 

120 

- 

140 

19 

170 

- 

175 

140 

- 

160 

20 

175 

- 

180 

160 

- 

180 

r 
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Since  the  results  of  the  analytic  solution  were  separable  in  all 
variables,  errors  were  also  separable. 

Consider  first  the  time  dependent  error. 


C(T) 


■ 

■-'‘j+i 

T(x)dx 

■T-Xj 

Y 

S(t)dt 

z 

j 

> 

L 

dx 

1 

S(t)T(T-t)dt 

100,  (10) 


Comparing  results  from  RENDER  and  the  analytic  solution  for  this  sample 
problem,  these  errors  were  found  to  be  negligible  before  times  around 
,1  msec  but  they  steadily  Increase  with  time  and  by  10  msec  have 
reached  82,86  per  cent.  This  error  is  shown  in  Fig,  11  which  is  the 
time  dependent  results  for  radial  Interval  between  125  to  150  meters  and 
deposition  angle  interval  between  80  to  90  degrees.  The  upper  curve  in 
the  figure  shows  data  from  RENDER,  the  lower  curve  shows  analytic  solution 
results.  To  verify  that  the  error  was  Indeed  a result  of  the  phenomenon 
described  by  Equation  10,  numerical  evaluation  was  made  with  the  sane 
time  bounds  {x^}  that  were  used  for  input  to  RENDER  and  a value  T • ,01 
second.  Evaluation  of  the  integral  over  the  source  term  only  in  Eq,  10 
was  made  at  the  end  points  since  the  source  was  analytically  Integrable; 


1 

S(t)dt  - ^ 
a 


-au. 


-au. 


(auj^+1)  e 


- (au2+l)  e 


(11) 


T-X 


J+1 


Dependent  Sample  Problem  Results  for  125^R^  150  Meters,  8O£0£  90  degrees 
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where  **2  " ^ ®^ror  value  of  84.20 

per  cent  resulted  from  this  evaluation  but  more  Important  was  the  fact 


that  for  the  first  49  elements  of  • the  evaluated  value  was  zero. 

Only  the  50th  element  of  t-Xj  (le.  .01  - .01  “ 0)  resulted  In  a non  zero 
term.  This  type  error  Is  Inherent  to  RENDER  and  will  be  manifested  when- 


ever the  larger  Intervals  of  transmission  time  develop  Intervals  In 
Eq.  11  that  contain  the  majority  of  the  source.  Fortunately,  those 
errors  that  are  Inherent  In  the  RENDER  results  which  are  angle  dependent 
will  typically  be  small  for  normal  production  problems. 


The  angle  dependencies  of  the  transmission  data  In  the  sample 
problem  were  such  that  very  large  errors  were  Introduced  into  the 
output  of  RENDER.  A comparison  of  some  Intermediate  RENDER  results  with 
the  limiting  function  of  the  analytic  solution  Is  shown  In  Fig.  12 
where  the  smooth  curve  Is  the  analytic  solution  and  the  histogram  Is 
from  Intermediate  data  In  RENDER.  This  normalized  histogram  Is  repre- 
sentative of  RENDER  results  for  all  radial  and  time  Intervals  when  the 
source  Is  between  0°  £ 90°.  As  shown  In  this  figure,  results  from 

RENDER  accurately  represent  the  shape  of  the  analytic  solution.  Careful 
checking  also  showed  the  magnitude  to  be  correct.  The  error  in  the 
RENDER  calculation  was  finally  Introduced  when  results  from  sources 
between  0 and  90  degrees  were  combined  with  results  from  sources  between 
90  and  180  degrees.  For  the  latter  set  of  data  the  transmission  angle 
set  Is  reversed.  Thus,  for  data  resulting  from  source  angles  between 
0 and  90  degrees  the  deposition  angle  set  was  {0,  20,  40,  60,  ...,  177, 
178,  179,  180}.  The  first  deposition  angle  Interval  In  the  final  printout 
of  RENDER  was  0 to  1 degree.  To  develop  data  for  that  Interval,  RENDER 
combined  the  data  from  sources  between  0 and  90  degrees  averaged  over  the 
Interval  0 to  1 degree  with  data  from  sources  between  90  and  180  degreea 
averaged  over  a 20  degree  Interval  (actually  0 to  20  degrees  In  final- 
result  space) . This  large-interval  small-interval  averaging  process  pro- 
duces errors  which  vary  directly  with  the  Interval  size  differential  and 
amount  of  variation  within  the  Intervals.  Comparisons  of  the  limiting 


ENERGY  (Normalized) 


DEPOSITION  POLAR  ANGLE  (deg.) 

Comparison  of  Total  Energy  Deposition  versus  Deposition  Polar 
Angle  as  given  by  RENDER  and  Analytic  Solution  for  Sample  Prob 
lem.  Source  Emission  for  Conical  Sources  between  0“  and  90® 


analytic  results,  the  analytic  results  averaged  over  deposition  polar 
angle  Interval  and  the  RENDER  results  are  shown  In  Fig.  13  for  the 
sample  problem  as  the  smooth  function,  the  histogram  without  circles 
and  the  histogram  with  circles  (when  it  differed  from  the  analytic 
histogram  results),  respectively.  It  is  seen  that  most  of  the  error 
between  the  RENDER  results  and  the  analytic  solution  occurs  at  deposi- 
tion polar  angles  near  0°  and  180°  where  Interpolation  on  the  transmis- 
sion data  was  required  In  RENDER. 

To  check  out  the  RENDER  program  using  the  energy  deposition 
data  computed  for  conical  sources,  a RENDER  problem  was  run  for  the 
gamma  ray  source  energy  interval  between  9 and  10  MeV.  The  gamma  ray 
source  was  taken  In  RENDER  to  be  Isotropic  with  a time  dependent  emission 
rate  described  by  the  source  function  given  in  Fig.  14.  The  deposition 
angle  bounds  and  time  values  used  to  print  the  RENDER  calculations  are 
given  in  Table  XIV. 

Table  XV.  lists  the  RENDER  results  integrated  over  time  as  a 
function  of  radial  position  and  deposition  angle  interval.  If  there  were 
no  error  in  the  RENDER  calculation,  then  the  energy  deposition  should  not 
vary  with  deposition  angle  interval  for  a given  radial  interval.  The 
data  in  Table  XV  indicate  that  there  were  some  errors  introduced  by 
RENDER  when  evaluating  the  energy  deposition  vs  radial  distance  for 
polar  angle  intervals  of  0 to  1 and  1 to  2 degrees.  For  polar  angle 
intervals  greater  than  2 degrees  there  is  very  little  difference  in  the 
energy  deposition  computed  for  each  deposition  polar  angle  interval  and 
each  radial  interval.  Also  shown  in  Table  XV  is  the  radial  distribution 
of  the  energy  deposition  for  a point  isotropic  source,  as  obtained  from 
the  DEPO  calculations  for  the  9-10  MeV  gamma  ray  conical  sources  that  were 
input  to  RENDER.  It  is  seen  that  the  RENDER  calculations  are  in  good  agree- 
ment with  the  DEPO  data  for  all  deposition  angle  intervals  greater  than 
2°  and  for  all  radial  intervals. 

The  time  dependent  energy  deposition  for  the  radial  interval  be- 
tween 330  and  400  meters  and  in  each  of  the  deposition  angle  intervals  is 


ENERGY  (Normalized) 


LIMITING  CASE  OF  ANALYTIC  RESULTS 


DISCRETE  ANALYTIC  RESULTS 


render  RESULTS 


DEPOSITION  POLAR  ANGLE  (Degrees) 


Comparison  of  Total  Energy  Deposition  given  by  RENDER  and  Analy' 
tic  Solution  for  Test  Problem:  Source  Emission  for  Conical 
Sources  between  0“  and  180“ 
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TABLE  XIV.  DEPOSITION  TIMES  AND  POLAR  ANGLE  BOUNDS  USED  IN  1 
RENDER  CALCULATION  FOR  A POINT  ISOTROPIC  9 TO  10 
CAMMA-RAY  SOURCE 

DEPOSITION  TIMES  (sec)  DEPOSITION  POLAR  ANGLE  BOUNDS 

1.000 

. 

10 

9.900 

9 

0 - 1 

5.000 

- 

10 

1.000 

- 

8 

1-2 

7.000 

- 

10 

1.000 

- 

8 

2-50 

8.000 

- 

10 

1.758 

- 

8 

50  - 90 

9.000 

- 

10 

2.330 

- 

8 

90  - 120 

1.000 

- 

9 

3.089 

- 

8 

120  - 150 

2.000 

- 

9 

4.095 

- 

8 

150  - 180 

3.000 

- 

9 

5.429 

- 

8 

3.500 

- 

9 

7.197 

- 

8 

4.000 

- 

9 

9.541 

- 

8 

4.500 

- 

9 

1.265 

- 

8 

5.000 

- 

9 

1.677 

- 

7 

5.500 

- 

9 

2.223 

- 

7 

6.000 

- 

9 

2.947 

- 

7 

6.  500 

- 

9 

3.907 

- 

7 

7.000 

- 

9 

5.180 

- 

7 

7.500 

- 

9 

6.867 

- 

7 

8.000 

- 

9 

9.102 

- 

7 

8.300 

- 

9 

1.207 

- 

6 

8.600 

- 

9 

1.600 

- 

6 

8.900 

- 

9 

2.121 

- 

6 

9.100 

- 

9 

2.812 

- 

6 

9.200 

- 

9 

3.728 

- 

6 

9.300 

- 

9 

4.942 

- 

6 

9.400 

- 

9 

6.551 

- 

6 

9.500 

- 

9 

8.685 

- 

6 

9.600 

- 

9 

1.151 

- 

5 

9.700 

- 

9 

1.526 

- 

5 

9.750 

- 

9 

2.024 

- 

5 

9.800 

- 

9 

2.683 

- 

5 

9.850 

- 

9 

3.557 

- 

5 
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— 

given  in  Table  XVI  for  times  less  than  5. 429X10  seconds.  Again  it  is 
seen  that  the  energy  deposition  rate  at  any  given  time  interval  for  de- 

r» 

position  polar  angle  intervals  less  than  2"  differs  somewhat  from  the  data 
for  deposition  polar  angle  intervals  greater  than  2°.  For  deposition 
polar  angles  greater  than  2°  there  is  very  little  variation  in  the  energy 
deposition  rate  in  a given  radial  Interval  with  deposition  polar  angle,  as 
expected . 

The  energy  deposited  in  the  atmosphere  for  print  times  less  than 

“8 

1X10  seconds  as  given  by  the  RENDER  calculations  was  found  to  vary  with 

time  in  the  same  manner  as  does  the  source  function  shown  in  Fig.  14. 

“8 

The  RENDER  results  for  times  less  than  10  seconds  can  thus  be  taken 

to  be  made  up  principally  of  the  energy  deposition  produced  by  first 

order  scattering  and  absorption  events.  The  results  of  integrating  the 

RENDER  data  for  150°  < [j/  ^ 180°  over  time  between  0 and  10  ^ sec  are 

given  in  Table  XVII  as  a function  of  the  radial  interval.  Also  shown  in 

Table  XVII  is  the  energy  deposition  data  for  an  isotropic  source  that  were 

obtained  from  the  conical  source  data  that  were  input  to  RENDER  for  9 to 

-20 

10  MeV  gamma  rays  in  the  0 to  10  sec  retarded  time  interval.  In  the 

conical  source  data  base  calculations,  the  energy  deposition  occurlng 
-20 

in  the  0 to  10  sec  retarded  time  interval  came  from  first  order  scat- 
tering and  absorption  events.  From  a comparlsion  of  the  data  in  Table 

XVII,  it  appears  that  there  is  a small  contribution  by  multiple  scattering 

—8 

to  the  RENDER  data  within  the  0 to  10  sec  retarded  time  Interval.  Hits 

contribution  results  from  the  fact  that  the  source  in  the  RENDER  calcu- 

"*8 

lation  emitted  radiation  over  a 10  sec  time  interval  whereas  in  the 
data  base  calculations  the  source  was  assumed  to  be  an  instantaneous 
emitter.  The  DEPO  time  dependent  energy  deposition  data  input  to  RENDER 
indicates  that  the  minimum  value  of  the  ratio  of  the  DEPO  to  RENDER  cal- 
culations shown  in  Table  XVII  should  not  exceed  0.91.  A review  of  the  data 
in  Table  XVII  reveals  that  this  ratio  is  always  greater  than  0.91,  as  expected. 

The  time  distribution  of  the  energy  deposition  in  the  600  to  800 

_g 

meter  radial  interval  for  times  greater  than  10  seconds  is  shown  in  Fig. 

15.  The  circled  points  in  Fig.  15  are  from  the  RENDER  calculation  and  the 
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TABLE  XVII.  COMPARISON  OF  RENDER  AND  DEPO  DATA  FOR  ENERGY 
DEPOSITION  RESULTING  FROM  FIRST  ORDER  SCATTFR- 
ING  AND  ABSORPTION  INTERACTIONS 


-3 

(keV  m /keV  of  source  energy') 

Radial  Interval 

(Meters) RENDER DEPO 


0 

10 

4.25-6* 

4.22-6 

10 

- 

20 

6.06-7 

6.04-7 

20 

- 

40 

1.39-7 

1.38-7 

40 

- 

60 

4.90-8 

4.84-8 

60 

- 

80 

2.51-8 

2.46-8 

80 

- 

100 

1.31-8 

1.28-8 

100 

- 

125 

8.61-9 

8.49-9 

125 

- 

150 

5.08-9 

4.99-9 

150 

- 

175 

3.57-9 

3.51-9 

175 

- 

200 

2.42-9 

2.53-9 

200 

- 

250 

1.53-9 

1.50-9 

250 

* 

300 

9.55-10 

9.31-10 

300 

- 

350 

6.08-10 

5.95-10 

350 

- 

400 

3.92-10 

3.80-10 

400 

- 

450 

2.62-10 

2.54-10 

4 50 

- 

500 

1.93-10 

1.87-10 

500 

- 

800 

1.20-10 

1.16-10 

600 

- 

800 

5.01-11 

4.82-11 

800 

- 

1000 

1.78-11 

1.71-11 

1000 

* 

1200 

6.76-12 

6.39-12 

1200 

- 

1500 

2.33-12 

2.18-12 

-6 


*Read  4.25-6  as  4.25X10 
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histogram  is  from  the  conical  source  data  input  to  RENDER.  There  is  good 

agreement  between  the  two  sets  of  energy  deposition  data  except  for  times 
— fl  —8 

of  /.iiXiO  sec  and  3.42X10  sec.  The  differences  at  these  times  is 

believed  to  result  from  the  fact  that  the  energy  deposition  data  input 

to  RENDER  was  for  an  instantaneous  emitting  source  and  the  source  used 

in  the  RENDER  calculation  emitted  over  the  time  interval  between  0 and 
“8 

10  sec.  Thus  the  convolution  in  RENDER  of  the  time  dependent  source 
data  with  the  time  dependent  energy  deposition  data  should  result  in  a 
slight  delay  of  the  energy  deposition  for  times  less  than  10  ^ sec. 

The  good  agreement  between  the  RENDER  calculations  and  the  conical 
source  data  input  to  RENDER  Indicates  the  RENDER  procedure  is  operating 
correctly.  The  results  of  the  test  problems  described  above  indicate  the 
possible  problems  that  might  arise  when  using  RENDER  to  convolute  a 
time,  energy  and  angle  dependent  source  term  with  the  nuclear  energy 
deposition  data  base  generated  for  conical  sources  in  this  study. 

The  nuclear  energy  deposition  data  base  and  the  RENDER  procedure 
are  operational  on  the  CDC-6600  computer  at  the  Los  Alamos  Scientific 
Laboratory  and  on  the  IBM-360/75  computer  at  the  Air  Force  Technical 
Applications  Center,  Patrick  Air  Force  Base,  FL.  The  RE2JDER  procedure 
was  written  in  FORTRAN-IV  language  for  each  of  these  computer  systems. 
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